Influence of magnetic on ferroelectric ordering in LuMnOs 
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We have studied the influence of antiferromagnetic ordering on the local dielectric moments of 
the MnOs and LuO? polyhedra by measuring neutron powder diffraction patterns of LuMnOa at 
temperatures near Tat. We show that the coupling is weak, because the magnetic exchange coupling 
is predominantly in the ah plane of the MnOs trigonal bipyramids, and the electric dipole moments, 
originating in the LuOy polyhedra, are oriented along the hexagonal c axis. Anomalies in the 
dielectric properties near Tn are thus caused by the geometric constraints between the MnOs and 
the LuO? polyhedra. 

PACS numbers: 



Multiferroic magnetoelectrics are materials that com- 
bine ferroelectric and magnetic ordering. These materials 
are of fundamental and technological interest as coupling 
between different order parameters enables multiple ways 
to interact with these materials, such as pressure, electric 
and magnetic fields. In most ionic materials the Coulomb 
repulsion between ions results in centred (i.e. non-polar) 
crystal structures. Conventionally, only for ions with lone 
pairs and non-magnetic ions, sizeable dipole moments 
can exist^. For few extraordinary materials magnetic 
and electric dipolar order can coexist. The coupling be- 
tween the magnetic and electric order must originate ei- 
ther from magnetic ordering affecting the electric dipoles 
via striction, or the displacements at the dipolar order- 
ing changing the magnetic order. These interactions de- 
pend intimately on the coupling between the magnetic 
and electric building blocks of the 3D-crystal structure. 

Recently the distorted perovskites BiMnO^Si^, 
BiFeO^^ and also ThMnOj^ have been shown to com- 
bine magnetic and ferroelectric order. Another class of 
materials that combine magnetic and electric dipolar 
order are the hexagonal manganites with general formula 
AMnOa, which are antiferromagnetic and ferroelectric. 
Hexagonal AMnOa (h-AMnOa), in particularly YMnOs, 
has attracted interest for applications, for instance in 
non-volatile memory devices® or as ferroelectric gate 
electrode^. Progress has also been made in the under- 
standing of the origin of the ferroelectric properties^i 
and the coupling between the ferroelectric and antiferro- 
magnetic orderi2ii£. Fiebig et al. have shown that below 
the antiferromagnetic ordering temperature T/v, the 
ferroelectric domain walls, or twin boundaries, always 
act as antiferromagnetic (AFM) domain walls as well. 
Furthermore, AFM domain walls also exist in the bulk 
of the ferroelectric (FE) domains. This suggests that not 
the orientation of the FE order and the AFM order are 
linked, but that the transition from one FE orientation 
to the other forces or accommodates a transition in 
the AFM ordesSiii. The coupling between electric and 
magnetic ordering requires detailed knowledge of the 
crystal structure. Unfortunately, there are conflicting 
results in the literature on the atomic positions when 



comparing neutron powder diffraction (NPD) and single 
crystal X-ray diffraction (SXD). Neutron diffraction re- 
sults show anomalous variations in the bond lengthsi^ii^ 
with temperature. 

The crystal structure of h-AMnOa has two major dif- 
ferences compared with the perovskite structure. Firstly, 
the Mn^^^ ions (high spin 3d'^) are located at the centre 
of a trigonal bipyramid instead of an octahedron. The 
resulting crystal field splits the 3c? levels in three energy 
levels, in order of increasing energy: xz and yz; xy and 
x^—y^; and iz^—r^. Consequently, Mn^+ 3c?^ has no par- 
tially filled degenerate levels and is therefore not Jahn- 
Teller activai^. Density of states calculations show some 
mixing between the Mn 3d and O 2p levelsi^ Secondly, 
the MnOs polyhedra are corner linked in sheets, sepa- 
rated by a layer of A ions. This contrasts the MnOg 3D 
network in perovskite manganites. As a result, h-AMnOa 
is pseudo-layered, reflected for instance in the Mn spin 
being in the xy plane even above T/v 

To explain the how the AFM and FE order are linked, 
we studied the temperature dependence of the crystal 
structure of LuMnOa around the antiferromagnetic or- 
dering temperature, T/v — 88 K, which is much lower 
than the FE ordering temperature, T > 573 Kii. Fur- 
thermore, we discuss the relation between the ionic radius 
of the A ion ta and the crystal structure. 



I. TEMPERATURE DEPENDENCE OF LuMnOg 

LuMnOa has been prepared with standard solid state 
synthesis. Neutron powder diffraction has been per- 
formed at the Polaris time-of-flight instrument at ISIS. 
The data have been refined using the GSAS package^?, 
including isotropic temperature factors. Sample quality 
has been checked by measuring a neutron powder diffrac- 
tion pattern at ambient temperature. No reflections from 
impurity phases could be detected and the refined struc- 
ture was in good agreement with single crystal data^^. 
Temperature dependent measurements, between 40 and 
120 K, were performed at the same beam line, using 
a standard Orange cryostat. Typically quality factors 
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for the refinements are = 2-963, R(F2) = 4.79% and 
Rt) = 1.74% for 47 variables. 



Fig. 1(a) shows the crystal structure of hexagonal man- 
ganite. It consists of non-connected layers of MnOs trig- 
onal bipyramids, that are corner-linked by in-plane oxy- 
gen ions, 03 and 04. The apical oxygen ions, 01 and 
02, form two close packed l ayers separated by a distorted 
layer of ions. In Fig. |l(b) a cross-section through 
the z ~ plane is given. The oxygen positions 01 and 
02 are similar to the Mn positions, albeit with z ~ ±i. 
The A site ions Al (A2) are located exactly above the 
positions 03 (04). In the next Mn ions containing layer 
{z ^ t he Mn, 01 and 02 positions are located at x 
(Fig. [1(b)]). 

The spacegroup of hexagonal LuMnOs in the ferroelec- 
tric phase is P63CTO. The hexagonal unit cell contains six 
formula units; there are seven inequivalent positions with 
10 refinable parameters. Since all symmetry elements are 
parallel to the c axis, all z coordinates are "free". For 
practical reasons, we have chosen to fix zmu =^ 0. The 
atomic positions are listed in Table 
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TABLE I: The atomic positions of hexagonal manganite 
AMnOa. m is the multiplicity of the site. 
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In Fig.|2]the lattice parameters has been plotted versus 
the temperature. The a axis decreases with decreasing 
temperature at a rate of 8 • 10^^ K~^. The lattice pa- 
rameter c is nearly independent of T, although a small 
discontinuity in c (~ 0.001 A) can be seen near T/y. SXD 
data at 295 K^^, a = 6.038(1) A, are in agreement with 
the extrapolation of Fig.|21 Thermal expansion data have 
been reported on powdered YMnOa samples at higher 
temperaturesSS. The expansion coefficient ~ 7 • 10~^ 
K"-"^ obtained from these data is one order of magnitude 
smaller than the value we find. 

Compared to the paraelectric centrosymmetric struc- 
ture, the MnOs bipyramids are rotated along an axis 
through the Mn ion and parallel to the 04-04 bond. This 
buckling moves the 03 (04) ions up (down) along the z- 
axis and moves the 01(02) ions at i, 0, ±^) along 
the X-axis. Fig. 13 shows the T dependence of the api- 
cal oxygen positions 01 and 02 of LuMnOa, plotted as 
the absolute difference with their paraelectric position at 
X = The X coordinate of these oxygen ions is a mea- 
sure for the buckling of the MnOs trigonal bipyramid, 
but small changes in xoi and 2:02 do not have a signifi- 
cant effect on the Mn-01 and Mn-02 bond lengths. The 
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FIG. 1: |(a)| View perpendicular to the c axis of the crystal 
structure of AMnOa. The trigonal bipyramids represent the 
MnOs and the circles the A i ons. This sketch highlights the 
stacking of the MnOs layers. |(b)| Sketch of the z=0 layer of 
LuMnOs, showing the two inequivalent Mn - O bonds. Mn, 
03 and 04 are shown as grey diamonds, closed and open cir- 
cles, respectively. Mn - 03 and Mn - 04 bonds are indicated 
by solid and dashed lines, respectively. The thin line outlines 
the unit cell. 

linear fits suggest that the buckling might reduce with 
decreasing T, but the slope and its error are the same 
order of magnitude. The angle of the line through 01 
and 02 with the xy-plane changes from 84.7° at T = 120 
K to 84.8° at T = 40 K. We conclude that there is hardly 
any change in the buckling of the MnOs layers with tem- 
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FIG. 2: Temperature dependence of lattice parameters of 
LuMnOs. 
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bond as is shown in Fig. 1(b) The in-plane Mn— 03 bond 
length is given by 



Mn-03 



-Mi 



(3) 



where in the approximation the effect of zqs 7^ can be 
ignored, since the partial derivatives are proportional to 
the magnitude of XMn and zos- The other in-plane bond 
length is given by 



Mn-04 = ^Ja^ixMn - 1/2)2 (^3/6)2} + c^z, 



2 

04- 



(4) 

Mn— 03 and Mn— 04 depend both on lattice parameter 
a and on the x coordinate of the Mn ion. We can write 
xm,i = I — (5, where S is the displacement of the Mn ion 
along the x axis from its paraelectric position at (i, 0, 0), 
see Fig. 0] 



Temperature (K) 




FIG. 3: Temperature dependence of the x coordinate of the 
apical oxygen positions of LuMnOa. The absolute difference 
with the paraelectric position at a: = i is plotted. Triangles 
and squares depict Ol and 02, respectively. The dashed lines 
are least squares linear fits. 
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The out-of-plane Al— 03 and A2— 04 bond lengths arc 
calculated using: 



Al-03 = c (zAi - Z03) 



(1) 



and 



FIG. 4: Sketch of the relation of the length of the two inequiv- 
alent in-plane Mn— O bonds, showing the effect of an in- plane 
Mn displacement. Symbols are the same as in Fig.|l(b)^ Thin 



dashed lines are lines parallel to the unit cell axes at 
tervals to guide the eye. 



A2-04 = c(zA2-^O4). (2) 

Since the temperature dependence of lattice parameter 
c is negligible, the T dependence of these bond lengths 
depend only on the T dependence of the z coordinates. 
We found in our refinements that the T dependencies 
of these z coordinates have strong positive correlations. 
Any changes in any of the z coordinates is mimicked 
closely by its companion. Therefore, the Lul— 03 and 
Lu2— 04 bond lengths show no dependence on tempera- 
ture. 



Then, applying xmu — ^ ~ Mn— 03 equals a(| — S) 
and Mn— 04 equals 

Mn-04 = y'a2{(l/6 + 5^ + 1/12} + c^z^^ - 

~ a (1/3 -I- 1/25), (5) 

where the effect of zqa 7^ and 5^ is negligible. 

From these approximate formulae we can see that a 
change, of —5a, in the length of the Mn— 03 bond length. 
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due to a displacement S of the Mn ion, is always accom- 
panied with a change in the length of the Mn— 04 bond 
of opposite sign and half the magnitude {+^5a). These 
in-plane Mn— O bond lengths are plotted against temper- 
ature in Fig. |5| Local extrema near Tpf are observed with 
the maximum value for Mn— 03 smaller than the mini- 
mum for Mn— 04. Below and above Tn the difference 
increases between the in-plane bond lengths. However 
the observed differences are not large, at most 0.041 A. 
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FIG. 5: Temperature dependence of in-plane Mn— O bond 
lengths around Tjv = 88 K. 

The out-of-plane Mn— O bond lengths are given by 



Mn-Oi 



xoiY + c^z^^ c zoi\ {i = 1, 2) 

(6) 

where the x components are neglected in the approx- 
imation, since xmu — xoi <SC zqi] (i = 1,2). The out-of- 
plane bond lengths also show local extrema near Tn, but 
the changes are less pronounced than for the in-plane 
bond lengths. The difference between the inequivalent 
out-of-plane Mn— O bond lengths is smaller than 0.02 A 
for aU T. 

The out-of-plane bond lengths are much shorter than 
the in-plane bond lengths, respectively ~ 1.87 A and 
~ 2.02 A. This is related to the crystal field splitting of 
the Mn 3d orbitals, with the unoccupied 3z^ — orbital 
in the out-of-plane direction. The antibonding character 
of the Mn 3d - O 2p interaction elongates the in-plane 
bond lengthsii . 

We note that with SXEli2i2iiS2i2^, it is found that the 
Mn— 04 bond length is smaller than the Mn— 03 bond 
length, whereas here we find that Mn— 04 is larger than 
Mn— 03. However, the difference between the inequiv- 
alent in-plane bond lengths is small (< 0.05 A). The 
out-of-plane Mn— O bond lengths are in line with the 
SXD data, and the difference between the inequivalent 
out-of-plane bond lengths is small (< 0.02 A) in agree- 
ment with the SXD data. Other NPD experiments on 
the crystal structure of hexagonal manganites in the lit- 
erature report larger differences between the in-plane or 



the out-of-plane Mn— O bond lengths. For instance, Kat- 
sufuji et al. report on LuMnOs at 300 K and find out- 
of-plane bond lengths of 1.98 A and 1.78 Ai^. Muiioz 
et al. find for YMnOa at room temperature a difference 
between in-plane bond lengths Ajn-pianc of 0.1 A^^. Also 
older structure reports using SXD show large differences 
between in-plane Mn— O bond lengths or out-of-plane 
Mn— O bond lengths. Yakel reports Ain-pianc = 0.08 A 
and Aout-of-piano = 0.09 A for LuMnOs^"*. Isobe re- 
ports Ai„_pia„c = 0.048 A and Aout-of-pianc = 0.058 A 
for YbMnOa at T 295 K^^. These SXD experi- 
ments however did not include reflections from all re- 
gions of hkl space, which is necessary for SXD on non- 
centrosymmetric crystals. 

In the previous section, we have discussed the tem- 
perature dependence of the bond lengths. Now, we will 
investigate whether the observed changes have any effect 
on the local dipole moments and the macroscopic ferro- 
electric moment. In the remainder of this article, we will 
define the local dipole moment associated with a cation 
site as the distance between the centre of gravity of the 
nearest neighbour O ions and the cation site. 

In Fig.|Slwe have shown that the in-plane Mn— O bond 
lengths show local extrema near Tjv, with the difference 
between them smallest at T/v . This means that also their 
contribution to the local dipole moments on the Mn site 
is smallest at T^. The space group symmetry relates 
the local dipole moment created by the small in-plane 
displacement of the Mn ion to the local dipole moments 
originating from the other two Mn ions on the same plane 
in the unit cell. Due to the sixfold screw axis, the local 
dipole moments are rotated over an angle of 120° with 
respect to each other. As a result, there cannot be any 
contribution to the macroscopic ferroelectric moment by 
the in-plane displacements. 

The difference between the out-of-plane Mn— 01 and 
Mn— 02 bond lengths create a small local dipole mo- 
ment < 0.02 A from the Mn coordination. These local 
dipole moments are parallel to the c axis and they are 
not cancelled by the crystal symmetry. However larger 
local dipole moments originate from the Lul and Lu2 
coordination**. The local dipole moments from the Lu 
coordination are at least an order of magnitude larger 
than the ones from the Mn coordination. 

The local dipole moments associated with the LuOs 
polyhedra are shown in Fig. The calculation of the 
centre of gravity of the oxygen polyhedra includes the 
seven ions at ~ 2.4 A and the one ion at ^ 3.2 A. Note 
that the local dipole moments from the Lu2 coordination 
are antiparallel with the ones from the Lul coordination. 

From Fig. it is clear that cooling through the Neel 
transition has no significant influence on the local dipole 
moments of LuMnOs. However, we have shown in Fig. [31 
that the local environment around the Mn ions does 
change. Since the Mn spins are in-planei^, even above 
T/v, the magnetic ordering transition at Tn is expected 
to have the largest effect on the in-plane Mn— O bond 
lengths. The ferroelectricity on the other hand depends 
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FIG. 6: Temperature dependence of LuMnOa of the local 
dipole moments from the Lul and the Lu2 polyhedra near 
Tat =88 K. Note that the right axis is reversed. 



FIG. 7: Lattice parameters for hexagonal AMnOa as a func- 
tion of TA (A=Lu, Er, Yb and Y) at T = 295 K. 



mostly on the out-of-plane Lu— O bond lengths^. Al- 
though the oxygen ions involved are the same, 03 and 
04, the Mn— O bond lengths change mostly by a dis- 
placement 8 of the Mn ion along the x axis, whereas the 
out-of-plane Lu— O bond lengths change because of an 
opposite movement of the Lul (Lu2) and 03(04) ions 
along the c axis. We conclude that there is no significant 
coupling between the ferroelectric moment of LuMnOa 
and the antiferromagnetic ordering temperature. This is 
in agreement with the conclusions from Ref. 0. 

II. RARE EARTH IONIC RADIUS 

In the second part of this article we will discuss the in- 
fluence of the rare earth ions, by their ionic radius r^, on 
the local dipole moments. The data used in this analysis 
has been taken from Refs. ^ (YMnOg); ^ (ErMnOg); 
m (YbMnOa) and ,19] (LuMnOg). 

In Fig. 13 the lattice parameters a and c a.t T = 295 
K are plotted against ta- The values of ta have been 
taken from Shannon and PrewittSSi. We observe that 
with increasing the lattice parameters increase. From 
LuMnOa to YMnOa lattice parameter a increases lin- 
early, with a total increase of ~ 1.7%. Lattice parameter 
c has a smaller increase of about 0.5%. We can clearly 
see from these data that the influence of ta is larger on 
the in-plane lattice parameter than on the out-of-plane 
lattice parameter. 

The crystal structure can be regarded as a stacking of 
hexagonally packed oxygen layers and A ion layers. One 
expects that 0—0 bond lengths are ^ 2.8 A, {tq — 1-4 
A) and A-0 bond lengths are ~ 2.4 A, {ta ^ 1.05 A). 
For instance, in LaMnOs the average 0-0 and La— O 
bond length (in the fcc-like closed packed 00 and AO 
layers) is ~ 2.7 A (r^a — 1-2 A). The observed values 
for the inter-plane A— 01 and A— 02 bond lengths in 
hexagonal AMnOs (not shown here) are consistent with 



the expected values. Likewise the observed increase with 
increasing rA of the A— 01 and A— 02 bond lengths (on 
average 1.6%) is in excellent agreement with the increase 
of a with rA- 

We note that in hexagonal AMnOa both the hexago- 
nally packed O layers at z « ±i and the MnO layers 
at z « have intra-plane 0-0 distances of ~ 3.6 A. 
Similarly, the in-plane A— A distances are also ^ 3.6 A. 
This is in contrast with the intra-plane bond lengths 
in pseudo-cubic perovskites. These anomalously large 
"bond lengths" can be understood, by considering the 
way the layers are stacked. Starting at the MnO layer 
at = 0, the stacking order is MnO - O - A - O - MnO 
with layer spacings of about -^c, j^c, j^c and ^c. One 
can visualise that the A ion layer (at z « i) is clamped 
between the closed packed O layer (at z « ^ and z w i). 
This will enlarge the spacing of the O and A ions in those 
layers. The larger the radius of the A ions, the more the O 
ions are pushed apart. However, the 0-0 bond length 
between the two neighbouring O layers at z « ^ and 
z ~ i is reduced to ~ 2.8 A, which corresponds to the 
value one expects (2 x tq). In this simple model, changes 
in rA can be accommodated completely by a change in 
the in-plane bond lengths and lattice parameter a. 

The decreasing distances in the ab plane with decreas- 
ing rA should also be reflected in the in-plane Mn— O 
bond lengths. The change of '--^ 1.6 % is in good agree- 
ment with the observed change of a with decreasing rA- 
Possible changes in the out-of-plane Mn-0 bonds are too 
small to be resolved within the scatter and error bars of 
the experiments. However, we observe a distinct relation 
between the Al— 03 and A2— 04 bond lengths and rA, 
as shown in Fig. |H1 

The short Al-03 and A2-04 bonds increase by 3.6% 
and 2.3% with increasing rA- Therefore, the "extra long" 
out-of-plane A— O bond lengths should decrease with in- 
creasing rA, since c is almost constant. Consequently, the 
dipoles of the A polyhedra decrease with increasing rA- 
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FIG. 8; out-of-plane A— O bond lengths as a function of ta- 



The increase in dipole moment from Y (ta ~ 1.075 A) 
to Lu {rA = 1.032 A) is about 10% as shown in Fig. 
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FIG. 9: Dipole moments of the Al and A2 polyhedra of h- 
AMnOs. Note that the A2 dipole is plotted on a reverse scale. 
Both dipoles tend to increase with decreasing rA- 

The transition from the paraelectric centrosymmetric 
to the ferroelectric non-centrosymmetric structure con- 
sists of two atomic displacements. The first is the buck- 
ling of the MnOs network, displacing 01 and 02 in the 
xy plane and 03 and 04 along the z axis. The second is 
the alternate displacement of the Al and A2 ions along 
the z axis away from the z = j mirror plane*. Formally, 
it has been shown that these atomic displacements con- 
sist of two modes. Firstly, there is a K3 mode, which 
changes the symmetry but can inherently not create a 
macroscopic ferroelectric moment. Secondly, there is a 
r mode, which does not change the space group symme- 
try, but introduces the macroscopic ferroelectricitjiSLSS. 
To see which of these two modes is infiuenced more by 
chemical substitution, we have plotted in Fig. ^| the z 



coordinates of the Al and the 03 sites. Both z coordi- 
nates of Al and 03 are about 0.02c away from their cen- 
trosymmetric equilibrium position. Clearly, the changes 
in the z coordinate of 03 site are much larger than the 
changes in the z coordinate of the Al site. The data in 
Fig. ^1 therefore show that upon chemical substitution 
of the A site the buckling increases with decreasing rA, 
which results in an increase of the deviation from the cen- 
trosymmetric positions. The displacement of the A ions 
from the mirror plane is more or less constant. This is 
supported by the changes in the x coordinates of the 01 
and 02 sites. Both are about 0.025a off their symmet- 
ric position at X = ^ and this deviation becomes larger 
by ^ 10% from R=Y to R=Lu. Wc conclude that the 
change in the local dipole moment is caused by atomic 
rearrangement and not by expansion of the unit cell. 
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FIG. 10; Fractional z coordinates of the Al and 03 position 
(0, 0, z) versus rA- Clearly the z coordinates deviate more 



from their equilibrium positions, respectively z 
0, with decreasing rA- 



and z 



The ferroelectric transition temperatures for these ma- 
terials have been reported, but different authors report 
completely different values or only "higher than" values. 
Furthermore, there are indications that the ferroelectric 
ordering temperature and the non-centrosymmetric to 
centrosymmetric phase transition temperature are dif- 
ferent by several hundreds degrees^S*^^. It is therefore 
not possible to relate the magnitude of the local dipole 
moments or the total dipole moment to the ferroelectric 
ordering temperature. 

The effect of changing can be summarised as fol- 
lows. A decrease in rA will reduce the "normal" A— O 
bond lengths. This will change the in-plane Mn— O bond 
lengths too. These changes are accompanied by an in- 
creased buckling of the MnOs bipyramids, parameterised 
by xoi, 2:02, zo3 and 204. The changes in the oxygen z 
coordinates induce an increase in the difference between 
the short and long Al— 03 and A2— 04 bond lengths. 
This increases the dipole moments from the Al and A2 
coordinations with decreasing rA- However, the antipar- 
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allel coupling between the local dipole moments at the 
Al and the A2 site prevents a clear effect on the total 
dipole moment. 

In conclusion, the hexagonal AMnOs compounds ex- 
hibit both antiferromagnetic and ferroelectric order. The 
magnetic exchange and spins are confined in the basal 
plane, whereas the dielectric moments originate predom- 
inantly from displacements parallel to the c-axis. Con- 



sequently, the distortions due to magnetostrictive effects 
and dipolar ordering are only weakly coupled. 
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